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Abstract 

The spectrum of neutrinos produced by the outflow of high temperature matter 
surrounding microscopic black holes is calculated for neutrino energies between one 
GeV and the Planck energy. The results may be applicable for the last few hours 
and minutes of a microscopic black hole’s lifetime. 
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1 Introduction 

Hawking radiation from black holes [jl| is of fundamental interest because it involves the 
interplay between quantum held theory and the strong held limit of general relativity. It is 
of astrophysical interest if black holes exist with sufficiently small mass that they explode 
rather than accrete matter and radiation. Since the Hawking temperature and black 
hole mass are related by = rrip/SrM, where mp = = 1.22 x 10^® GeV is the 

Planck mass with natural units h = c = /cB = l,a present-day black hole will evaporate 
and eventually explode only if Tu > 2.7 K (microwave background temperature). This 
requires M < 4.6 x 10^^ kg which is approximately 1% of the mass of the Earth. More 
massive black holes are cooler and therefore will absorb more matter and radiation than 
they emit, hence grow with time. Taking into account emission of gravitons, photons, 
and neutrinos a critical mass black hole today has a Schwarszchild radius of 68 microns 
and a lifetime of 2 x 10^^ years. 
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There is uncertainty expressed in the literature about whether the radiation emitted 
by such small black holes is emitted independently, or whether there is sufficient inter¬ 
action among the emitted quanta to form a partially thermalized and rapidly expanding 
fluid around the black hole. Heckler has given arguments in favor of the latter for Hawking 
temperatures above 1 GeV |@]. Cline, Mostoslavsky and Servant solved the Boltzmann 
equation in the relaxation-time approximation for such temperatures and found that sig- 
nihcant particle scattering did occur, although not enough for perfect fluid flow. One 
of us showed that a self-consistent description of an outgoing fluid just marginally kept 
in local equilibrium could be given, but it required the assumption of sufficient initial 
particle interaction, and viscosity played a crucial role [^]. That paper was followed by an 
extensive numerical analysis of the relativistic viscous fluid equations by us [§]. We also 
calculated the spectrum of high energy gamma rays expected during the last days, hours, 
and minutes of the black hole’s lifetime. We suggested that the most promising route for 
discovery of such microscopic black holes is to search for point sources emitting gamma 
rays of ever-increasing energy until suddenly the source shuts off. 

In this paper, a follow up to our previous work, we focus on high energy neutrino 
emission from black holes with Hawking temperatures greater than 100 GeV and corre¬ 
sponding masses less than 10® kg. It is at these and higher temperatures that new physics 
will arise. Such a study is especially important in the context of high energy neutrino 
detectors under construction or planned for the future. Previous notable studies in this 
area have been carried out by MacGibbon and Webber 0 and by Halzen, Keszthelyi and 
Zas 0 , who calculated the instantaneous and time-integrated spectra of neutrinos arising 
from the decay of quark and gluon jets. 

The source of neutrinos in the viscous fluid picture is quite varied. Neutrinos should 
stay in thermal equilibrium, along with all other elementary particles, when the local 
temperature is above 100 GeV. The reason is that at energies above the electroweak scale 
of 100 GeV neutrinos should have interaction cross sections similar to those of all other 
particles. Thus the neutrino-sphere, where the neutrinos decouple, ought to exist where 
the local temperature falls below 100 GeV. The spectra of these direct neutrinos are 
calculated in section H. 

Pions and muons remain in local thermal equilibrium down to temperatures on the 
order of 100 to 140 MeV, as shown in |^. Neutrinos also come from decays involving these 
particles. The relevant processes are (i) a thermal pion decays into a muon and muon- 
neutrino, followed by the muon decay /x —>■ ez/eix^, and (ii) a thermal muon decays in the 
same way. The spectra of neutrinos arising from pion decay are calculated in section HI 
while those arising from direct or indirect muon decay are calculated in section IV. 

The spectra from all of these sources are compared graphically in section V. We also 
compare with the spectra of neutrinos emitted directly as Hawking radiation without 
any subsequent interactions. The main result is that the time-integrated direct Hawking 
spectrum falls at high energy as E~^ whereas the time-integrated neutrino spectrum 
coming from a fluid or from pion and muon decays all fall as E~^. Thus the fluid picture 
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predicts more neutrinos at lower energies than the direct Hawking emission picture. If a 
microscopic black hole is near enough the instantaneous spectrum could be measured, and 
its shape and magnitude would provide information on the number of degrees of freedom 
in the nature on mass scales exceeding 100 GeV. Conclusions are presented in section VI. 

Any details not provided here on the outflow of high temperature matter surrounding 
a black hole may be found in 


2 Directly Emitted Neutrinos 

In this section we hrst review the emission of neutrinos by the Hawking mechanism un- 
modihed by any rescattering. Then we estimate the spectra of neutrinos which rescatter 
in the hot matter. The last scattering surface should be represented approximately by 
that radius where the temperature has dropped to 100 GeV. The reason is that neutrinos 
with energies much higher than that have elastic and inelastic cross sections that are com¬ 
parable to the cross sections of quarks, gluons, electrons, muons, and tan leptons. Much 
below that energy the relevant cross sections are greatly suppressed by the mass of the 
exchanged vector bosons, the W and Z. Furthermore, the electroweak symmetry is broken 
below temperatures of this order, making it natural to place the last scattering surface 
there. A much better treatment would require the solution of transport equations for the 
neutrinos, an effort that is perhaps not yet justihed. All the formulas in this section refer 
to one flavor of neutrino. Our current understanding is that there are three flavors, each 
available as a particle or antiparticle. The sum total of all neutrinos would then be a 
factor of 6 larger than the formulas presented here. 


2.1 Direct neutrinos 


The emission of neutrinos by the Hawking mechanism is usually calculated on the basis 
of detailed balance. It involves a thermal flux of neutrinos incident on a black hole. 
The Dirac equation is solved and the absorption coefficient is computed. This involves 


numerical calculations 
by 


rm . The number emitted per unit time per unit energy is given 


T, 1 

dEdt 271 expi^E/Tn) + 1 ’ 

where is an energy-dependent absorption coefficient. There is no simple analytic 
formula for it. For our purposes it is sufficient to parametrize the numerical results. A 
fair representation is given by 


r. 


27 ^2 / 0.925 

Tjj y exp(5 — aE/Tn) + 1 


( 2 ) 
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and a ~ 1.607. The exact expression has very small amplitude oscillations arising from 
the essentially black disk character of the black hole. The parametrized form does not 
have these oscillations, but otherwise is accurate to within about 5%. 

If there are no new degrees of freedom present in nature, other than those already 
known, then the time dependence of the black hole mass and temperature are easily found. 
The temperature is 


ThW 


1 /— 

Stt \ Sat ) 


(3) 


Here the temperature is Tnito) = Tq at (negative) time to and increases to inhnity at 
time t = 0. The constant a is approximately 0.0043 for Th > 100 GeV. This relationship 
between the time and the temperature allows us to compute the time-integrated spectrum, 
starting from the moment when Th = Tq. There is a one-dimensional integral to be done 
numerically. 


dNf'^ 27ml 0.925 \ 

dE 16{47r)^aE^ Jo exp(a;) + 1 \ exp(5 — ax) + 1/ 


(4) 


In the high energy limit, meaning E ^ Tq, the upper limit can be taken to inhnity with 
the result 


29.9mp mp 
dE ^ {ATT)^aE^ ~ 565^' 


2.2 Direct neutrinos from an expanding fluid 

Neutrinos emitted from the decoupling surface have a Fermi distribution in the local rest 
frame of the huid. The phase space density is 

qE'/T^ ^ I ■ 

The decoupling temperature of neutrinos is denoted by T^. The energy appearing here is 
related to the energy as measured in the rest frame of the black hole and to the angle of 
emission relative to the radial vector by 


E' = 7 y(l — v,y cos 9)E . 


(7) 


No neutrinos will emerge if the angle is greater than 7r/2. Therefore the instantaneous 
distribution is 


dEdt ‘'(2irV 


d(cos 6) cos 0f\E, cos 6) 





T 
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exp[-nE{'y^ - u^)/Tl\ + 


T 
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( 8 ) 
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where is the radius of the decoupling suface and u = v'y. Integration over the energy 
gives the luminosity (per neutrino). 

We need to know how the radius and radial flow velocity at neutrino decoupling 
depend on the Hawking temperature or, equivalently, the black hole mass; we have already 
argued that the neutrino temperature at decoupling is about = 100 GeV. Our numerical 
solutions to the viscous fluid equations ^ resulted in a simple scaling law between the 
flow velocity and radius. 


u{r) = Us 



(9) 


Here rs is the Schwarszchild radius and us ~ 0.10. (The viscous fluid results may not 
be applicable until the radius exceeds rs several times over.) Hence is known if is. 
Specihcally 


1 

4:7iTh \us) 


( 10 ) 


The hnal piece of information is to recognize that each type of neutrino will contribute 
7/8 effective bosonic degree of freedom to the total number of effective bosonic degrees of 
freedom in the energy density of the fluid. After integrating the instantaneous neutrino 
energy spectrum to obtain its luminosity we equate it with the appropriate fraction of the 
total luminosity. 


^fluid _ 7 3-V^ ^ ^Hawking 

8 90 106.75 


( 11 ) 


where 


^Hawking _ 647r2«T^ . (12) 

The number 106.75 counts all effective bosonic degrees of freedom in the standard model 
excepting gravitons. Without the graviton contribution a is approximately 0.0042. This 
results in an equation which determines v, equivalently u, in terms of the Hawking tem¬ 
perature. 


(3 - Vu)vtul 
(1 - 


2i3457r^^ g 


427 


-auc 




(13) 


Numerically the constant C = 1.2 x 10“®. 

We are interested in black holes with = 100 GeV. The corresponding range 

of neutrino-sphere flow velocities corresponds to 0.2 < < oo. There is no simple 

analytical expression for the solution to eq. (13) for this wide range of the variable. At 
asymptotically high temperatures the left side has the limit 16u^. For intermediate values 
the left side may be approximated by 22u1. We approximate the left side by the former 
for Uu > 11/8 and by the latter for 0.2 < < 11/8. Thus 




22T4 ) 


for Ui, < 11/8 


(14) 
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1/8 


for Ui, > 11/8 . 


(15) 


(CJt\ 

Ii6tJ 

This approximation is valid to better than 20% within the range mentioned. 

The time-integrated spectrum can be calculated on the basis of eqs. (8), (3), (10) and 
either (13), which is exact, or with the approximation of (14-15). The resulting formulas 
are lengthy and not very interesting, therefore not displayed here. What is interesting is 
the high energy limit, E /:§> lOTj, = 1 TeV. This is determined wholly by the hrst of the 
two exponentials in eq. (8) with the coefficient of 1. Using also eqs. (3), (10) and (14) 
we hnd the limit 

dN^nid 4185C(6)C'1/^ m|,T^ ^ 

' 85U% ^ ^ 1040U4 ^ 

This E~^ spectrum is characteristic of all neutrino sources in the viscous fluid description 
of the microscopic black hole wind. 


3 Neutrinos from Pion Decay 

The decoupling temperature Tf for pions was estimated in to be in the range from 
100 to 140 MeV, comparable to the pion mass and substantially less than the neutrino 
decoupling temperature. Muon-type neutrinos will come from the decays of charged pions, 
namely —>■ and tt~ —>■ In what follows we calculate the spectrum of muon 

neutrinos; the spectrum for muon anti-neutrinos is of course identical. 

In the rest frame of the decaying pion the muon and the neutrino have momentum q 
determined by energy and momentum conservation. 





(17) 


Numerically q = 0.2134m^. 
malized to one, is 


In the pion rest frame the spectrum for the neutrino, nor- 


6(E' - q) 
(Ep' Anq 


(18) 


The Lorentz invariant rate of emission is obtained by folding together the spectrum with 
the rate of emission of tt’*'. 


(Epdt 



1 ^ f EE^ - P • P' 

47rg V mjr 


(19) 


The spectrum of pions is computed in the same way as the spectrum of direct neutrinos 
from the expanding fluid in sect. 2.2 but with one difference and one simplification. The 
difference is that the pion has a Bose distribution as opposed to the Fermi distribution 
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of the neutrino. The simplihcation is that the fluid at pion decoupling has a highly 
relativistic flow velocity with Uf ^ 'jf 1. The analog to eq. (8) is 



( 20 ) 


The instantaneous energy spectrum of the neutrino is reduced to a single integral. 



( 21 ) 


Here T'min = + q^)/2Eq is the minimum pion energy that will produce a neutrino 

with energy E. We are interested only in high energy neutrinos with E S> in which 
case E yr^\r^ = m^T'/2q' is an excellent approximation. The integral can also be expressed 
as an inhnite summation 



( 22 ) 


We have found from the numerical solutions that 



(23) 


and that 



(24) 


Here df = 12 is the effective number of bosonic degrees of freedom at the decoupling 
temperature Tj. These results take into account the total luminosity of the black hole 
minus that contributed by gravitons and neutrinos. When E S> JjTf only the hrst term 
in the summation in eq. (22) is important. 

The instantaneous spectrum is easily integrated over time because the flow velocity at 
pion decoupling is highly relativistic. The spectrum arising from the last moments when 
the Hawking temperature exceeds Tq is 



dxx^ e (25) 


In the limit that E S> 7 /(To)Tj the upper limit on the integral may be taken to inhnity 
with the following result. 


diV; 4647r 203^ ^^^ m|>T/ ^ 



(26) 


dE ^ 427 \m^) \ 2df ) E^ ^ 300E^ 


This is much smaller than the direct neutrino emission from the fluid because the neutrino 
decoupling temperature is much greater than the pion decoupling temperature Tf. 
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4 Neutrinos from Muon Decay 


Muons can be emitted directly or indirectly by the weak decay of pions: tt” —/r“z/eh^ 
plus the charge conjugated decay. Both sources contribute to the neutrino spectrum. The 
invariant distribution of the electron neutrino (to be specific) in the rest frame of the 


muon is 



4 

Timt 


(3m^ 


AE')E' 


(27) 


where the electron mass has been neglected in comparison to the muon mass. This 
distribution is used in place of the delta-function distribution of eq. (18) for both direct 
and indirect muons, as evaluated in the following two subsections. 


4.1 Neutrinos from direct muons 

The instantaneous spectrum of Ve, t'e, or arising from muons in thermal equilibrium 
until the decoupling temperature of T/ can be computed by folding together the spectrum 
of muons together with the decay spectrum of neutrinos in the same way as eq. (19) was 
obtained. Using eqs. (20) and (27) results in 


dEdt 


37 ^ 7 / \ 


nE 

9 -^ + 2 

27 / 1 / 



-|- exp 



' 107 / 1 / 

nE 


-7-2 


nE 

T/U. 



(28) 


where Ei is the exponential-integral function. In the high energy limit, dehned here by 
E 3> 27 /T/, the spectrum simplihes to 


dWy- 2r}T,E / _ 

dEdt Sn'jf y 2'^fTf) 


(29) 


The time-integrated spectrum can be calculated in a fashion analogous to that followed 
in section 3. Thus 

dNA ^ _ 870 / 203(7 \ mpdf 

dE Qln^df \ 2df ) E^ 
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In the high energy limit this simplihes to 


11745 

dE 4887r®(i( 


[225C(7) - 217C(6)] 


_ mlTf 

2df / ^ 250^4 


4.2 Neutrinos from indirect muons 

The spectrum of neutrinos coming from the decay of muons which themselves came from 
the decay of pious proceeds exactly as in the previous subsection, but with the replacement 
of the direct muon energy spectrum by the indirect muon energy spectrum. In the rest 
frame of the pion the muon distribution is 


E'—^ 

d^p' 


/ 


6{p' - q) , 


which is the hnite mass version of eq. (18). Folding together this spectrum with the 
spectrum of pions (20) yields the spectrum of indirect muons. 


dEdt 


m^rjTf “1 ( 

- 2^ ^ exp - 

'Trn ' ^ 1 


nrriTrE 


This spectrum is now folded with the decay distribution of neutrinos to obtain 


dEdt 


m^rjTf “ 1 r / 27 2 2 3 ■ 

- -—- > — < Ei[—nx) 5- n x - n x' 

V 6 3 

/ 19 23 2 2 2\ -nxV^ 

-1- nx + -rcx^ e > 

V 6 6 3 /1 


where 


a;± = 


E r(m^ + ± g 


27 /T/ 

The high energy limit E ^ '^fTf is 


^2^mdirM m^rjT] fe-^- e-^+1 

dEdt 37rg x1_ x\ 

The integral over time can be done in the usual way. The resulting expression cannot 
be expressed in closed form, is lengthy, and is not illuminating. The high energy limit is 
simple and of the familiar form 1/E^. 

^jymdirM ^74^ / 203(2 ^+ q^ (m^ + 2g2) 


2989df \ 2df 


1250.^4 
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5 Comparison of Neutrino Sources 

In this section we compare the different sources of neutrinos that were computed in the 
previous sections. All the hgures presented display one type of neutrino or anti-neutrino. 
That type should be clear from the context. For example, equal numbers of z/g, t'g, 
z/^, Vti and Vr are produced as Hawking radiation and by direct emission by the fluid 
at the neutrino decoupling temperature T^. Only electron and muon type neutrinos are 
produced by muon decay, and only muon type neutrinos by pion decay. These differences 
could help to distinguish the decay of a microscopic black hole from other neutrino sources 
if one happened to be within a detectable distance. 

The instantaneous spectra are displayed in Fig. 1 for a Hawking temperature of 1 
TeV corresponding to a black hole mass of 10^ kg and a lifetime of 7.7 minutes. The 
instantaneous spectra for a Hawking temperature of 10 TeV corresponding to a black 
hole mass of 10® kg and a lifetime of 0.5 seconds are displayed in Fig. 2. There are 
several important features of these spectra. One feature is that the spectrum of direct 
neutrinos emitted by the fluid, at the decoupling temperature of = 100 GeV, peaks at 
a lower energy than the spectrum of neutrinos that would be emitted directly as Hawking 
radiation. The peaks are located approximately at y/T^Tn for the fluid and at 5Th for 
the Hawking neutrinos. The reason is that the viscous flow degrades the average energy 
of particles composing the fluid, but the number of particles is greater as a consequence 
of energy conservation. In the viscous fluid picture of the black hole explosion direct 
neutrinos are emitted as Hawking radiation without any rescattering when Th < T^, 
whereas when Th > Ty they are assumed to rescatter and then be emitted from the 
neutrino-sphere located at Ty. It is incorrect to add the two curves shown in these hgures. 
In reality, of course, it would be better to use neutrino transport equations to describe 
what happens when Th ~ Ty. That is well beyond the scope of this paper, and perhaps 
worth doing only if and when there is some observational evidence for microscopic black 
holes. 

Another feature of Figs. 1 and 2 to note is that the average energy of neutrinos arising 
from pion and muon decay is much less than that of directly emitted neutrinos. Again, 
the culprit is viscous huid how degrading the average energy of particles, here the pion 
and muon, until the time of their decoupling at T/ ^ 100 MeV. On the other hand their 
number is greatly increased on account of energy conservation. The average energies of 
the neutrinos are somewhat less than their parent pions and muons because energy must 
be shared among the decay products. The spectrum of muon-neutrinos coming from the 
decay tt —pz/ is the softest because the pion and muon masses are very close, leaving 
very little energy for the neutrino. 

The time-integrated spectra, starting at the moment when the Hawking temperature 
is 1 and 10 TeV, are shown in Figs. 3 and 4, respectively. The relative magnitudes and 
average energies reflect the trends seen in Figs. 1 and 2. At high energy the Hawking 
spectrum is proportional to E~^ while all the others are proportional to E~*, as was 
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already pointed out in the previous sections. Obviously the greatest number of neutrinos 
by far are emitted at energies less than 100 GeV. The basic reason is that only about 5% 
of the total luminosity of the black hole is emitted directly as neutrinos. About 32% goes 
into neutrinos coming from pion and muon decay, about 24% goes into photons, with 
most of the remainder going into electrons and positrons. 


6 Observability of the Neutrino Flux 


We now turn to the possibility of observing neutrinos from a microscopic black hole 
directly. Obviously this depends on a number of factors, such as the distance to the black 
hole, the size of the neutrino detector, the efficiency of detecting neutrinos as a function 
of neutrino type and energy, how long the detector looks at the black hole before it is 
gone, and so on. 

For the sake of discussion, let us assume that one is interested in neutrinos with energy 
greater than 10 GeV and that the observational time is the last 7.7 minutes of the black 
hole’s existence when its Hawking temperature is 1 TeV and above. As may be seen 
from Fig. 3, most of the neutrinos will come from the decay of directly emitted muons. 
Integration of eq. (31) from i^min = 10 GeV to inhnity, and multiplying by 4 to account 
for both electron and muon type neutrinos and anti-neutrinos, results in the total number 
of 


N. 


m%Tf oi 

^ ^ 8 X 10^^. 


750 


(38) 


This does not take into account neutrinos directly emitted from the fluid. For i^min = 1 
TeV, for example, eq. (16) should be used in place of eq. (31) (see Fig. 3), and taking 
into account tan-type neutrinos too then yields a total number of about 3 x 10^®. For an 
exploding black hole located a distance d from Earth the number of neutrinos per unit 
area is 


N^{E > 10 GeV) 
K{E>1 TeV) 



km 


-2 



km ^ . 


(39) 

(40) 


Although the latter luminosity is smaller by three orders of magnitude, it has two advan¬ 
tages. First, 1/3 of that luminosity comes from tan-type neutrinos. Unlike electron and 
muon-type neutrinos, the tan-type is not produced by the decays of pions produced by 
interactions of high energy cosmic rays with matter or with the microwave background 
radiation. Hence it would seem to be a much more characteristic signal of exploding black 
holes than any other cosmic source (assuming no oscillations between the tan-type and the 
other two species). Second, the tan-type neutrinos come from near the neutrino-sphere, 
thus probing physics at a temperature of order 100 GeV much more directly than the 
other types of neutrinos. 
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What is the local rate density piocai of exploding black holes? This is, of course, 
unknown since no one has ever knowingly observed a black hole explosion. The hrst 
observational limit was determined by Page and Hawking |^. They found that the local 
rate density is less than 1 to 10 per cubic parsec per year on the basis of diffuse gamma 
rays with energies on the order of 100 MeV. This limit has not been lowered very much 
during the intervening twenty-hve years. For example, Wright used EGRET data to 
search for an anisotropic high-lattitude component of diffuse gamma rays in the energy 
range from 30 MeV to 100 GeV as a signal for steady emission of microscopic black holes. 
He concluded that piocai is less than about 0.4 per cubic parsec per year. If the actual rate 
density is anything close to these upper limits the frequency of a high energy neutrino 
detector seeing a black hole explosion ought to be around one per year. 


7 Conclusion 

This paper has been a continuation of our previous work on a viscous fluid description of 
the radiation from microscopic black holes. In this paper we have calculated the spectra 
of all three flavors of neutrinos arising from direct emission from the fluid at the neutrino- 
sphere and from the decay of pions and muons from their decoupling at much larger 
radii and smaller temperatures. We should emphasize in particular the usefulness of 
distinguishing between the electron and muon-type neutrinos and the tan-type ones. The 
latter are much less likely to be produced by high energy cosmic rays. If the rate density 
of exploding microscopic black holes in our vicinity is anywhere close to the current limit 
based on gamma rays, it should be possible to observe them with present and planned 
large astrophysical neutrino detectors. Observation of high energy neutrinos, especially in 
conjunction with high energy gamma rays, may provide a window on physics well beyond 
the TeV scale. 
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Figure 1: The instantaneous direct neutrino spectra emerging from the fluid with neu¬ 
trino-sphere located at Ty = 100 GeV compared to the direct Hawking radiation. Also 
shown are neutrinos arising from pion, direct muon and indirect muon decays at a decou¬ 
pling temperature oiTf = 100 MeV. Here the black hole temperature is Th = 1 TeV. All 
curves are for one flavor of neutrino or anti-neutrino. 
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Figure 2: Same as Fig. 1 but with = 1 TeV. 
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Figure 3: The time integrated neutrino spectra emerging from a microscopic black hole. 
Here the calculation begins when the black hole temperature is Tq = 1 TeV. Either the di¬ 
rect Hawking radiation of neutrinos or the direct neutrino emission from a neutrino-sphere 
at a temperature of 100 GeV should be used. All curves are for one species of neutrino 
or anti-neutrino. 
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Figure 4: Same as Fig. 3 but with Tq = 10 TeV. 
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